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Abstract 

We report a comprehensive density functional theory (DFT) + U study of the energetics of 
charged and neutral oxygen defects in both PuC>2 and a-Pu203, and present a quantitative deter- 
mination of the equilibrium compositions of reduced Pu02 (PuC>2— x ) as functions of environmen- 
tal temperature and partial pressure of oxygen, which shows fairly agreement with corresponding 
high-temperature experiments. Under ambient conditions, the endothermic reduction of Pu02 to 
CK-PU2O3 is found to be facilitated by accompanying volume expansion of Pu02-x and the possi- 
ble migration of O-vacancy, whereas further reduction of a-Pu203 is predicted to be much more 
difficult. In contrast to the endothermic oxidation of PuC>2, the oxidation of a-Pu2C>3 is a stable 
exothermic process. 
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I. INTRODUCTION 



Plutonium oxides have long been of great interest not only for their high importance in 
nuclear fuel cycle and in long-term storage of Pu-based radioactive waste, but also their 
fancy physical properties and complex chemical reactions. Most of the rich behaviors of 
Pu ions in different oxidation states can be attributed to the complex character of Pu-5/ 



orbitals 



U-y], which locates in the boundary of localized and delocalized 5/ states among the 
actinides. Since different oxidation states have widely different solubilities and demonstrate 
distinct effects on chemical corrosion of plutonium surface, the redox properties of Pu-oxides 
play a crucial role in the plutonium science and the reactor safety. 

When exposed to oxygen-rich air, metallic Pu surface oxidizes readily to a protective 
dioxide (PUO2) layer{4-6], with a thin layer of the sesquioxide (PU2O3) existing in between. 
Although, under special aqueous condition, the hydration reaction on Pu0 2 surfaces (PuOo-f- 
xH 2 0— 7-Pu02+a; + xH 2 ) has been reported to generate higher compound Pu0 2+x (x<0.27)|7|, 
it has been theoretically proved to be strongly endothermicjs] and the products (Pu0 2+:c ) are 
experimentally found to be chemically unstable at elevated temperatures {jj. Thus, Pu0 2 
is the highest stable oxide as the compound of choice for the long-term storage of Pu. 
However, under oxygen- lean conditions (in the vacuum or inert gas), the Pu0 2 -layer can 
be reduced to sesquioxides (Pu 2 3 ), which can promote the hydrogenation of Pu metal 4, 
[lo| . Besides the well-known sesquioxide /3-Pu 2 C>3 in hexagonal structure (P3ml), the low- 
temperature sesquioxide phase is cubic a-Pu 2 C>3 (ia3), with ideal stoichiometric structure 
similar to the 2x2x2 fluorite Pu0 2 (Fm3m) supercell containing 25% O vacancy located in 
the 16c (0.25,0.25,0.25) sites [see Fig. l(top)]. In fact the cubic sesquioxide was detected 
in an abnormal body centered cubic PuOi.5 2 form 11], since a-Pu 2 03 is stable only below 
600 K and usually, a mixture of cubic PuOi 52 and PuOi. 98 can be obtained by partial 



reduction of Pu0 2 at high temperature and cooling to room temperature 



12|. Meanwhile, 



the complex Pu-0 phase diagram[13j shows that the Pu0 2 -Pu 2 C>3 region mainly consists 
of widely nonstoichiometric compounds Pu0 2 _ x and has not been fully described since the 
nonstoichiometry can remarkably vary the properties of Pu oxides at high temperatures. 
As the essential component of the uranium-plutonium mixed oxide (MOX) fuels in fast re- 



actors, considerable efforts have been dedicated to the defect chemistry of Pu oxides 



and theoretical models have been proposed based on the measured properties/parameters 
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such as the formation energy of oxygen defects. More recently, Minamoto et al.\\j\ 
employed the DFT calculated E* of single neutral O-vacancy to discuss the phase equi- 
librium of PUO2-X-PU2O3, however, the conventional DFT that applies the local density 
approximation (LDA) or generalized gradient approximation (GGA) underestimates the 
strong correlation of the 5/ electrons and, consequently, describes PuO"2 as incorrect fer- 
romagnetic FM conductor instead of antiferromagnetic AFM Mott insulator reported 



by experiment 
U 



20]]. Fortunately, several beyond-DFT approaches, the DFT (LDA/GGA) + 



n 



211 ] . the self-interaction corrected LDA|2|, |22|, the hybrid density functional of (Heyd, 



Scuseria, and Enzerhof) HSE[23|, and LDA + dynamical mean-field theory (DMFT) 24 1 
have been developed to correct the failures of conventional DFT in calculations of Pu oxides. 
The effective modification of pur e DFT by LDA/GGA+f/ formalisms has been confirmed 



27-. 



35|, which have reported the AFM- insulator ground state of 



widely in studies of Pu02 
Pu0 2 . 

In this work, we present a comprehensive DFT+£7 study of the charge states and for- 
mation energies of oxygen point defects in both Pu02 and a-Pu203, and the equilibrium 
compositions of Pu02- X as functions of the environmental parameters (temperature and 
partial pressure of oxygen). 

The rest of this paper is organized as follows. The details of our calculations are described 
in Sec. 2. In Sec. 3, we present and discuss the results. In Sec. 4, we summarize our main 
conclusions. 



II. DETAILS OF CALCULATION 



DFT+U studies 
Dudarev et al. 
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The calculations are carried out using the Vienna ab initio simulation package |25| in the 
framework of the projector-augmented wave (PAW) method 26]. Following our previous 
27H32J of plutonium oxides, we apply the GGA+U scheme formulated by 
to correct the strong on-site Coulomb repulsion among the Pu-5/ elec- 
trons. The choice of U and J parameters are based on an overall agreement between available 
experimental data and theoretical results as regarding the basic physical properties of PUO2 
and PU2O3, including the lattice constants, bulk modulus, phonon spectra and density of 
states (DOS), and the valence electronic DOS. Based on our extensive test calculations, it 



is found that the Hubbard U can alter the electronic-structure properties of Pu02 as well as 
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those of CK-PU2O3. Specifically, when the effective parameter U e //(=C/-J) is less than 2.0 eV, 
the calculated electronic DOS indicates that a-Pu203 are metallic FM-conductor instead of 
the Mott-insulators. And the combination of £7=4.75 and J=0.75 eV is proved to allow the 
ground-state properties of both PuC>2 and a-Pu203 to be reasonably described. A 2x2x2 
supercell derived from the ideal fluorite PuC>2 is used, which accommodates 64 O and 32 Pu 
with the collinear lk AFM order. Following the recent DFT+t/ calculations of the ground 



state of lk AFM insulating U0 2 



361 ]. we have carefully tested the use of such collinear lk 



AFM order to make the convergence to the ground state systematically, and strictly abided 
by the basic principle that the state with lowest total energy corresponds to the ground- 
state. The charged oxygen defects are modeled by adding or removing electrons from the 



CSS 



supercell. The geometry optimization is carried out until the forces on all atoms are 
than 0.01 eV/A. The plane -wave cut-off energy of 400 eV and a 7x7x7 Monkhorst-Pack 
/c-point mesh (64 irreducible k points) are applied, and the convergences with respect to 
these input parameters are checked to be within 0.001 eV. 

The formation energy of an isolated O-vacancy El or interstitial E{ in charge state q can 
be defined as 



E s v = E vac - E per + n (T,po 2 ) + qE F , (1) 

E( = E int - E per - hq(T, p 02 ) + qEp, (2) 

where E vac , E int , and E per represent the total energies of the system with an O-vacancy, 
an O-interstitial, and the perfect system, respectively. Ep is the Fermi energy relative to the 
valence band maximum (VBM). fio{T,po 2 ) is t ne oxygen chemical potential under partial 
pressure po 2 . Since the available experiments have just provided the related data [in Fig. 
4(b)] under the conditions of such high temperature and particularly low pressure, where the 
interaction between the gaseous molecules is negligibly small, the ideal gas approximation of 
gaseous Oojs expected to be reasonable. Thus, we can use the well-known thermodynamic 



expression [381] 



Ho(T,po 2 ) = - {Eo 2 +fi 02 (T,p°) + k B T\n(po 2 /p )) , (3) 
where Eq 2 is the total energy of the oxygen molecule. For the standard pressure p° — 1 



atm, the values of jlo 2 (T,p ) have been tabulated in Ref. 39[. With the u ab initio atomic 
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FIG. 1: (Color online) Top: crystal structures of PuO"2 and a-Pu203 and their building blocks with 
grey arrows indicate the directions of the atomic displacements. The blue, red and white spheres 
denote the Pu, O atoms and O vacancies. Bottom: the atom-projected density of states (PDOSs) 
for Pu-5/ and 0-2p electrons in PuC-2 (left) and a-Pu203 (right). The Fermi level is set to eV. 



thermodynamics" method 381]. one can see that the zero-temperature DFT+U calculations 
just provide input parameters (total energies) for the thermodynamic formalism of the redox 
energetics [equations (1) and (2)], and the environmental dependence of El and E{ is there- 



fore dominated by fio, which includes three required and correlated factors, namely, T, po 



and no 2 . When compared to the experimental value of 2.56 eV 40j, the DFT-GGA ca 



2 



cula- 



tion is known to overestimate the binding energy of molecular O2 by 0.8 eV/0.5O2(27|, |41| . 
and the correction term of 0.8 eV is thus applied in this work. The calculated oxidation re- 
action energies in Pu 2 03+|0 2 — > Pu0 2 are —3.81 and —3.19 eV for /3-Pu 2 3 and a-Pu 2 3 , 
respectively, which indicate that a-Pu203 is more stable at K. 



III. RESULTS AND DISCUSSION 



We begin by an overall comparison of the atomic- and electronic-structure properties 
between Pu02 and a-Pu203 presented in Fig. 1. In Pu02, each Pu cation is bonded to 
eight O anions and every O is embedded in a regular Pu-cation tetrahedron, whereas after 
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an isostructure reduction to a-Pu203, the coordination of Pu-cation reduces to six and O 












) 
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nm is close to the experimental values ranging from 1.103 to 1.107 nmjll 
a little larger than 2<2q u02 of 1.093 nm (in theory) and 1.079 nm (in experiment 7j), then 
an interesting volume expansion (7.6% in theory vs 7.4% in experiment) appears during 
the PuC>2 — » a-Pu2C>3 isostructure reduction. Moreover, in PuC^-x, the abnormal volume 
expansion is found to mount up almost linearly with the increasing nonstoichiometry x 
[see Fig. 4(a) below] in good agreement with experiments. The orbital-resolved atom- 
projected density of states (PDOS) suggests that the reason for such an abnormal issue is 
the localization of Pu-5/ electrons after the formation of O- vacancy. Specifically, the highest 
occupied band (HOB) with a range of —5 to eV is mainly the 5/(Pu)-2p(0) hybridization 
in PuC>2. Whereas in a-Pu203, the width of HOB decreases to 1.2 eV and the Pu-5/ and 
0-2p occupied bands are well separated. The deduced populations of 5/ state are 4.1 and 5.3 
electrons in Pu02 and CK-PU2O3, respectively. Thus, the formation of O-vacancy in a-Pu203 
or reduced Pu02 withdraws 5/ electrons from the bonding state and causes the coulomb 
repulsions among Pu cations with more localized 5/ electrons. Both Pu02 and a-Pu203 are 
insulators with the calculated gaps of 1.8 (the same as experiment [6j ) and 2.0 eV. 

For the isolated oxygen defects of different charge states in Pu02 and a-Pu203, the 
calculated zero-temperature El and E\ at the O-rich limit (i.e. po 2 =p°=l atm) as functions 
of the Fermi energy Ep in the insulating band gaps are shown in Fig. 2. Depending on 
the location of Ep, the thermodynamically stable charge states of O vacancies are found 
to vary between two values of +2 and in Pu02 [see Fig. 2(a)], and transit from +2 to 
via +1 in a-Pu 2 3 [Fig. 2(c)]. In both Pu0 2 and a-Pu 2 3 , it appears that the lowest El 
corresponds most of the time to the charge state of +2. And according to the signs and 
values of El, one can conclude that the reduction of Pu0 2 (to a-Pu 2 3 ) is endothermic, 
and the further reduction of a-Pu203 will be much more difficult due to the much higher 
values of El. During the reduction of Pu02 to nonstoichiometric PUO2-2; compounds, the 
defect-chemistry study by Stan et al. [16] predicted that +2 O-vacancy is the dominating 
defect species at T = 1373 K in the very low nonstoichiometry (x < 0.01) region, namely, 
+2 O-vacancy dominates only at the beginning of the redox process, then +1 O-vacancy in 
a limited intermediate region, and finally neutral O-vacancy in the deep region (x > 0.025). 
In current calculations based on a 2x2x2 Pu02 supercell, the nonstoichiometry x is not 
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FIG. 2: (Color online) Formation energy E-F of oxygen point defects in different charge states as a 
function of Fermi energy in the band gap of Pu02 and a-Pu203. 

less than 0.03125 so that neutral O-vacancy will be the only dominating defect species 
according to the foregoing prediction. In both oxides, O interstitials can be —2 or charged 
depending on the position of Ep. In contrast to the endothermic oxidation of Pu02 (to 
higher oxides Pu02+ X ), a minus E{ indicates that the oxidation of a-Pu203 (to PUO2) is a 
stable exothermic process. 

Assuming that both Pu-oxides are in equilibrium with an external environment and 
translating /io into the environmental conditions (T and po 2 ) according to Eq. (3), the 
results presented in Fig. 3(a)-(d) mainly show how a specific ideal gas approximation corrects 
the formation energies and influences the critical oxygen partial pressure p c Q2 with E£=0 or 
E{ =0 eV, namely, the evolutions of formation energies of neutral O defects as a function of 
T and po 2 ■ The O point defects are found to be sensibly modulated by the environmental 
T and po 2 as follows: (i) both E? and E{ are independent of po 2 when T = 0K; (ii) the El 
linearly decreases with decreasing ln(po 2 ) when T > K and reaches zero at a critical Pq 2 , 
which increases with increasing temperature; (iii) the E\ linearly decreases with increasing 
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FIG. 3: (Color online) Dependence of the calculated formation energy of oxygen point defects in 
Pu02 and a-Pu203 on temperature T and oxygen partial pressures po 2 . 

ln(po 2 ) when T > K and as T is increased, the critical p c 02 decreases for Pu0 2 , however 
increases due to the minus value of E{ for a-Pu 2 3 ; (iv) when ln(p O2 )=0 (i.e., po 2 =p Q =\ 
atm), under the ideal-gas approximation the fio 2 (T ^ 0,p°) term can notably depress the 
but increase E( when compared with the corresponding formation energy at T = K, and 
it is particularly so at high temperatures. Thus, the fio 2 (T ^ 0, p°) term is a contributory 
factor to figure out the critical p c 02 values, since it includes contributions from rotations and 
vibrations of molecular O2, as well as the ideal-gas entropy at 1 atm. 

Now we focus on the environment- dependent equilibrium compositions of nonstoichio- 
metric Pu0 2 _x- The four points of p c Q2 (at T=300, 400, 600, and 900 K) in Fig. 
3(a), correspond to four environmental conditions, under which the PuOi. 96875 com- 
pound is thermodynamically stable. For the further discussion of other PuC^-z com- 
pounds, the averaged formation energy El of neutral O-vacancy can be written as 

— No l^vac — E per + N _vfi (T,po 2 )} , where N _y is the total number of O-vacancy in 
the Pu32 64-7v o _ v supercell (with x= N °~ v ). According to the atomic structure of a-Pu203 
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FIG. 4: (Color online) (a) The average formation energy E v of O-vacancy (left axis) and the 
optimized volume V of PuG^-a; (right axis) varies as the x. Solid and open squares - &l with and 
without volume optimization; Blue circles - optimized V; Blue spheres - experimental V (compiled 



from Refs 
Refs. 



42] and 



3) 



45j). (b) The calculated (spheres) and experimental (squares, compiled from 



4514481]) relationships of equilibrium composition of PuG^-z vs temperature and oxygen 



partial pressure. Lines are used to guide the eyes. 

(see Fig. 1), we first adopt the Pu32064-at o _ v supercell with even (iVo-v = 2, 4, 6, 8, 12, 
and 16) vacancies located along the parallel (111) diagonals. In Fig. 4(a), the calculated 
&f (with and without volume optimization) and the optimized volume V of Pu 32 64 _iVo_v; 
together with the compiled experimental V of PUO2-X, are plotted as functions of x. Al- 
though the volume expansion is found to effectively depress the E~l when x > 0.25, the El 
appears to monotonically increase with x and when x=0.5 it reaches 4.71 eV, which is much 
larger than the reduction energy E r of 3.19 eV in Pu02 =>■ CV-PU2O3 + |02- This is due to 
the special and fully optimized distribution of O-vacancy in the stable a-Pu2C>3, specifically, 
vacancies located along four unequal and nonintersecting (111) diagonals. We then partially 
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optimize the O-vacancy distribution based on the Pu3 2 064-tv o _ v supercell with optimized 
V (for iVo-v= 4, 6, 8, 12), and the corresponding El is found to decrease with the increas- 
ing x and tend to approach E r of 3.19 eV. To accomplish the optimum distribution, there 
may be many complex and correlated mechanisms and processes, among which the possible 
migration of O-vacancy and the accompanying volume expansion of the PuO^-x matrix are 
doubtless very important and will be studied in our next work. 

Figure 4(b) shows equilibrium compositions of Pu0 2 _ x reported by this work as functions 
of environmental T and po 2 , which appear to be in good agreement with experimental data 
at high temperatures of 1173 K, 1373 K and 1773 K. Thus, under special high-temperature 
and low-pressure conditions, the ideal-gas approximation is a reasonable enough request for 
a quantitative determination of the environment (T and po 2 ) dependent phase equilibrium 
of Pu0 2 _ x . As a comparison, since the flo 2 (T,p ) term under the ideal-gas approxima- 



tion is not considered, the equi 
recent theoretical studies 



43 



ibrium compositions of Ce0 2 _ x and Hf0 2 ± x reported by 
44j are found to be only qualitative agreement (instead of 
quantitative) with the corresponding high-temperature experiments. At low-temperature 
conditions, our results (not shown here) indicate that it is very hard to achieve the desired 
thermo-equilibrium conditions. However, with a reducing atmosphere such as CO or H 2 , 
these conditions can be generated. Furthermore, in our DFT+[7 calculations, it is found 
that some cluster structures of O vacancies can efficiently decrease the average formation en- 
ergy E£, the fact of which indicates that at low temperatures the optimized configuration of 
O vacancies in the Pu0 2 _ x matrix is the actual key factor that influences the corresponding 
thermo-equilibrium conditions, whereas at high temperatures the O vacancies potentially 
have an even distribution. 



IV. CONCLUSIONS 

In summary, based on the DFT+[7 framework and the ideal-gas approximation, we sys- 
tematically investigate the environmental dependence of the redox energetics of Pu0 2 and 
a-Pu 2 03. Our results clearly reveal that the reduction of Pu0 2 is an endothermic process, 
whereas the reduction of a-Pu 2 03 will be much more difficult due to the higher formation 
energies of O-vacancy. In contrast to the endothermic oxidation of Pu0 2 , a minus formation 
energy of O-interstitial indicates that the oxidation of a-Pu 2 03 is exothermic. And the re- 
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ported equilibrium compositions of PuC^-x show quantitative agreement with experiments. 
Our current study may provide a guiding line to the understanding of redox behaviors of 
plutonium oxides. 
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